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I report on a study of finite temperature QCD by the CP-PACS Collaboration toward 
a precise determination of the equation of state with dynamical u,d quarks. Based on a 
systematic simulation using improved Wilson-type quarks on lattices with temporal size 
Nf = 4 and 6, the energy density and pressure are calculated as functions of temperature 
and renormalized light quark mass in the range T/T^ ~ 0.7-2.5 and mps/my = 0.65- 

0. 95. Results for Nt = 4 are found to contain significant scaling violations, while results 
for Nt = 6 are suggested to be not far from the continuum limit. On the other hand, the 
quark mass dependence in the EOS turned out to be small for mps/my ;^0.8. 

1. INTRODUCTION 

Theoretical understanding on the nature of the finite temperature QCD phase transition 
and its thermodynamic properties is indispensable in extracting an unambiguous signal 
of quark-gluon plasma production from heavy-ion collision experiments. In particular, 
the equation of state (EOS) of quark-gluon plasma belongs to the most basic information 
needed. Here, numerical simulations of lattice QCD provides us with the only systematic 
way to calculate these quantities directly from the first principles of QCD |]1|. 

In this paper, I present our study of finite temperature QCD on the lattice with dynami- 
cal u and d quarks, focusing on the efforts toward a precise determination of the EOS [0,0. 
Calculations are made on the CP-PACS computer, a dedicated parallel computer devel- 
oped at the University of Tsukuba in 1996 by physicists and computer scientists [|]. With 
2048 node processors the peak performance of the CP-PACS is 614.4 GFLOPS. Intensive 
lattice QCD calculations on the CP-PACS have clarified the existence of quenched errors 
in the traditional calculations of light hadron spectrum in the quenched approximation , 
and opened the way towards precise full QCD calculations of hadronic quantities through 
the first systematic study of the light hadron spectrum with dynamical u,d quarks [Q. 
Our simulations of finite temperature QCD are based on the experiences of these zero- 
temperature studies. 

2. FINITE TEMPERATURE QCD ON THE LATTICE 
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Figure 1. (a) Phase diagram and simulation points on 16^ x 4, 16^ x 6 and 16^ lattices, (b) 
Lines of constant physics and of constant temperature. Solid lines are mps/my constant 
lines, and dashed lines are T /Tp^ constant lines for A^^ = 4. The values of T /Tp^ for the 
dashed lines are given on the right edge of the figure. 



2.1. Lattice fermions 

Until recently, EOS with dynamical quarks has been computed only with the Kogut- 
Susskind (staggered) type lattice fermions because of a smaller computational demand 
and persistence of a part of the chiral symmetry on the lattice |^ . The staggered fermion 
action, however, only allows four degenerate quark flavors, and the physically interesting 
cases of 2 and 2 + 1 flavors are simulated by artificially modifying a coefficient in the 
update algorithm, effectively introducing a non-local action. It has also been found that 
the critical scaling for two-flavor chiral transition extracted with this formalism does not 
reproduce the theoretically expected 0(4) critical exponents 

These problems of the staggered fermions make it imperative to study the issue with 
alternative lattice fermions — Wilson-type fermions or recently proposed lattice chiral 
fermions. The latter fermions, however, require too much computational cost on current 
computers to perform a systematic investigation over a wide range of the parameter space. 

Wilson-type fermions have the manifest flavor symmetry and locality for any number of 
flavors. On the other hand, chiral symmetry is explicitly broken at finite lattice spacings, 
which complicates the phase diagram analysis on one hand and introduce sizable lattice 
artifacts on the other hand. Here, improvement of the lattice action was shown to be 
effective in reducing the lattice artifacts in finite temperature QCD [j^. In particular, the 
expected 0(4) scaling was observed around the two-flavor chiral transition |1T2[ . 



From these observations, and also from a preparatory study of improved actions 
we adopt an RG improved gauge action coupled with a clover-improved Wilson quark 
action, as in our zero-temperature full QCD simulations [Q. 

2.2. Phase diagram and lines of constant physics 

Our simulation parameters are summarized in |@J^. We study A^j = 4 and 6 lattices. 
The temperature is given by T = 1/Nta, where we determine the lattice spacing a by 
the vector meson mass my at T = 0. Results for the phase diagram, together with the 
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simulation points for the EOS study, are summarized in Fig. [T|(a). The horizontal axis 
is the inverse bare gauge coupling (3 = Q/g"^. The vertical axis, the hopping parameter 
K, corresponds to the freedom of the bare quark mass. The chiral limit is given by the 
Kc{T = 0) line. The triangular region to the left, surrounded by the lines Kc{Nt = 4), 



is the parity-flavor broken phase []10| , pnj| for = 4. Dashed lines Kf are the finite 
temperature pseudocritical lines for A^^ = 4 (left) and 6 (right), which are determined by 
the peak position of the susceptibility for a Z(3)-rotated Polyakov loop. 

In previous studies of EOS with staggered-type quark actions, the pressure and energy 
density are determined as functions of temperature for a fixed value of bare quark mass 



m 



^"■^'^a or m'i°'^^/T = m^f'^'^Nta. While m'i°-^'^a and Nt are practically easy to set in simula- 



tions, the results at different temperatures represent values for different physical systems 
with different renormalized quark masses. This is not useful for phenomenological appli- 
cations; we want to see the temperature dependence of thermodynamic observables for 
a fixed physical system with fixed renormalized quark mass, i.e., on a line of constant 
physics. 

Therefore, it is important to evaluate the lines of constant physics also in finite tem- 
perature physics. Our results for the lines of constant physics, which are defined by fixed 
values of the ratio mps/my of the pseudo-scalar to vector meson masses at zero tempera- 
ture, are summarized in Fig. |l|(b) by solid lines. In the same figure, we also plot the lines 
of constant temperature normalized by the pseudocritical temperature Tpc at Kt{Nt = 4) 
on the same line of constant physics. 

3. RESULTS 

3.1. Chiral transition temperature 

We first study the chiral transition temperature in the limit of massless quarks. In 
the phase diagram. Fig. |I](a), the chiral transition point is the crossing point of the chiral 
limit line Kc{T = 0) and the finite temperature pseudocritical line Kt. We confirmed that 
the 0(4) critical scaling, expected from an effective sigma model, is well satisfied with 
our data 0. This implies that the chiral transition is second order in two-flavor QCD. 

From a theoretical argument, we expect that the parity-broken phase at flnite Nt locates 



in the low temperature phase [11]. This provides us with an estimate of the lower bound 
for the chiral transition point; jSct > 1.538(46) from the Kc{Nt = 4) lines in Fig. |T](a). We 
can also estimate (3ct using the 0(4) scaling properties of observables: From an 0(4) flt to 
the chiral condensate, we flnd (3^ = 1.469(73). On the other hand, from an 0(4) scaling 
ansatz to the Kt line, we obtain P^t = 1.557(28). Collecting all the results, and converting 
to the physical units using the vector meson mass at T = 0, we obtain = 171 ± 4 MeV 
as our best estimate 0]. 

3.2. Equation of state 



We calculate the pressure p by the integral method |TJ]. For the energy density e, we 
combine the results for p and those for the interaction measure e — 3p obtained by the 
differential method. Our results are shown in Fig. |^. 

The continuum limit corresponds to the limit of large A^^. Our data show a 50% decrease 
from A^i = 4 to 6, both in the pressure and energy density, which is too large to attempt 
a continuum extrapolation. On Nt = 6 lattices, however, the magnitude and temperature 
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• mp5/m„=0.65 
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► mpjm„=0.95 



Q- 4.0 




2.5 3.0 



25.0 



10.0 



5.0 



0.5 




SB N|=6 



SB continuum 

• mpg/m„=0.65 
■ mpj./m„=0.70 

♦ mpj./m„=0.75 
▲ mp5/m„=0.80 

mp3/m„=0.85 
T mpg/m„=0.90 
► mp_/m„=0.95 



1.0 



1.5 



T/T„, 



2.0 



2.5 



3.0 



Figure 2. Pressure (a) and energy density (b) on 16^ x 4 (filled symbols) and 16^ x 6 
(open symbols) lattices, as functions of T/Tpc [||. 



dependence of the two quantities are quite similar between our improved Wilson quark 



action and the staggered quark action [|T^. Together with the fact that the Nt = Q results 
are close to the continuum SB limit at high temperatures, the approximate agreement of 
EOS between two different types of actions may be suggesting that the Nt = 6 results are 
not far from the continuum limit. 

In Fig. 0, different values of mps/my correspond to different renormalized light quark 
masses. We find that the dependence on the quark mass is quite small for mps/my ~ 0.65- 
0.8. A weak quark mass dependence appears only at mps/my ^0.9 for A^^ = 4 (the errors 
for the Nf = 6 data are still large to conclude a quark mass dependence). We find that 
the pressure and energy density in the quenched QCD [mps/my = 1) are much smaller 
[0; only about 1/7 of the values at mps/my ~ 0.7. See Figs. 15 and 17 of [|]. These 
results suggest that most dynamical quark effects are saturated between mps/my = 1 
and 0.9. Therefore, although the light u,d quark mass point, mps/my = m.„/mp = 1.8, 
is still far away, our values for the EOS at mps/my ~ 0.65-0.8 may be already close to 
those at the physical point, except in the very vicinity of the chiral transition point where 
a singular limit according to the 0(4) critical exponents is expected. 

This result of small quark mass dependence may not be surprising since hadron mass 
results in our zero-temperature simulations show that the renormalized quark mass 
at /i = 2 GeV in the MS scheme at mp^/my ^ 0.65-0.8 equals m^'^{^ = 2GeV) 
100-200 MeV, which is similar in magnitude to the critical temperature Tc ~ 170 MeV 
estimated for two-flavor QCD. For comparison, finite mass corrections for free fermion 
gas only amount to 7% when the temperature equals the fermion mass rrif, and exceed 
50% only when mj/T^3. 

4. CONCLUSIONS 

We have shown that a study of finite temperature QCD with Wilson-type quarks is 
feasible when improved lattice actions are adopted. Based on a systematic determination 
of the lines of constant physics, we calculate thermodynamic quantities in improved lattice 
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QCD with two flavors of Wilson-type quarks on iVt = 4 and 6 lattices, in the range 
n^'Ps/'iT^v = 0.65-0.95 and T/Tc ^2.5. We found that the quark mass dependence is small 
for mps/mv ^0.8. 

Concerning the continuum extrapolation, we find large scaling violations in EOS at 
Nf = 4, while Nt = 6 results are suggested to be not far from the continuum limit. 
Therefore, a precise continuum extrapolation may be possible from simulations at A^^ ^6. 
This is, however, computationally quite demanding. In this connection, I would like to 
introduce our recent study of EOS using anisotropic lattices [0: We found that, in 
quenched QCD, anisotropic lattice is efficient to reduce discretization errors in EOS, such 
that a precise continuum extrapolation is possible with much smaller amount of computer 
time. Anisotropic lattices may help performing continuum extrapolations also in QCD 
with dynamical quarks. 

The studies presented in this paper are performed by the CP-PACS Collaboration. 
This work is in part supported by the Grant-in- Aid of Ministry of Education, Science and 
Culture (Nos. 12304011 and 13640260 ). 
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